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Abstract. In this short note we provide two extensions on the recent explicit 
results on the matrix-product ansatz for the non-equilibrium steady state of a 
markovianly boundary-driven anisotropic Heisenberg XXZ spin 1/2 chain. We write a 
perturbative solution for the steady state density matrix in the system-batyh coupling 
for an arbitrary (asymmetric) set of four spin-flip rates at the two chain ends, 
generalizing the symmetric-driving ansatz of [Phys. Rev. Lett. 106, 217206 (2011)]. 
Furthermore, we generalize the exact (non-perturbative) form of the steady state 
for just two Lindblad channels (spin-up flipping on the left, and spin-down 
flipping on the right) to an arbitrary (asymmetric) ratio of the spin flipping rates 
[Phys. Rev. Lett. 107, 137201 (2011)]. In addition, we also indicate a simple proof of 
uniqueness of our steady states. 



PACS numbers: 02.30.Ik, 03.65.Yz, 05.60.Gg, 75.10.Pq 
1. Introduction 

The anisotropic Heisenberg (XXZ) model ^ of n coupled quantum spins 1/2 with the 
Hamiltonian 

^ = E (2-;^i + 2^<i + ^^H+i) (1) 

i=i 

can be considered as a prototype of a many body quantum model of strong interactions. 
We write Pauli operators on a tensor product space Tn = (C^)®", as (j| = ® o"^ (g) 
l2"-j, Id being a d-dimensional unit matrix, where cr^ = ^{a"^ ± icr^) and a'"'^'^ are the 
standard 2x2 Pauli matrices. 

The XXZ model exhibits a rich variety of equilibrium and nonequilibrium physical 
behaviors. In nature it provides an excellent description of the so called spin-chain 
materials [2], and it is believed to provide the key for understanding of various 
collective quantum phenomena in low dimensional strongly interacting systems, such 
as magnetic or superconducting transitions in two dimensions. Although equilibrium 
(thermodynamic) properties of XXZ chain are well understood in terms of Bethe Ansatz 
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[5] , as the model represents a paradigmatic example of quantum integrable systems, its 
nonequilibrium properties at finite temperature are a subject to a lively debate 

We consider the non-equilibrium quantum transport model [5] based on a 
Markovian master equation for the XXZ chain in the Lindblad form [6l [7] 

^ = -i[H, p{t)] + 2L,p{t)Ll - {LlL,, Pit)} (2) 

k 

with a set of four boundary-supported Lindblad operators 

Li = y/acr^, L2 = \fp>ot, = y/icr~, L4 = v^a-+ (3) 

where a,/3,j,6, represent, respectively, spin-down/spin- up incoherent transition rates 
at the left and right boundary of the chain. These are assumed to be the only incoherent 
processes in the model, whereas its bulk dynamics is fully specified by the Hermitian 
many-body Hamiltonian H. 

In particular we are interested in the non- equilibrium steady state (NESS), with the 



density operator poo = hmt^oo pit) satisfying the fixed point condition: 

Cpoo = 0, (4) 

where C is the Liouvillean, decomposed into the unitary, and dissipative parts: 

C:= -iadH + 1), (5) 

i) : = aV^- + + 7^.- + SV^^ , (6) 

Aip) :=2apat-{ata,p}. (7) 



The boundary driven open XXZ model ([2]) with ffTllSl) can be considered as a 
hybrid classical/quantum Markovian lattice gas (see e.g. [8j for related ideas), namely 
the boundary injection/absorbtion rates are the same as in some classical Markovian 
stochastic many-body processes (i.e. simple exclusion processes [9l [10]) while the 
bulk dynamics is fully coherent. Note that the effect of incoherent processes on the 
boundary, which works against developing strong macroscopic entanglement in the 
course of time-evolution, also enables efficient applicability of Liouville space density- 
matrix-renormalization group methods for computing NESSes for generic (non-solvable) 
local spin chain Hamiltonians [11]. The applicability of incoherent boundary processes 
to model the (magnetic) baths can be indeed justified if there is a finite correlation 
(coherence) length in the microscopic model of the baths. 

Below we make a few remarks concerning recently developed exact solution for 
NESS of the boundary driven XXZ model [T31[T5]. In section [2] we show how uniqueness 
of NESS simply follows from a theorem of Evans [12] • In section [3] we then extend the 
zeroth and the first order of NESS density matrix in the weak-coupling perturbative 
expansion [13] to the case of arbitrary boundary spin-fiipping rates a,/3,7,5, while in 
section m we in a similar way extend the non-perturbative exact solution of Ref. [15]. 
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2. Proof of uniqueness of the steady state 

Let us first show that under quite general conditions, the open- XXZ model (E]) with 
f|T|3|) . possesses a unique NESS ([5]), i.e. the fixed point poo is independent of the initial 
state p(0). 

We start by noting a theorem of Evans [12] (see also |13|) which essentially 
states that NESS is unique iff the set of operators {iJ, Li, L2, . . .} generates, under 
multiplication and addition, the entire algebra of (bounded) operators, in our case the 
Pauli algebra B{J-'n) of the spin-1/2 chain on n sites. Indeed, this is true even if we take 
only the Hamiltonian H, and a single pair of one up-flip and one down-flip Lindblad 
operators out of four ([3]), say af, and crjf. Note that the scalar prefactors 1/0,1/^ 
are not important for this discussion as we are only interested in the generators of the 
algebras and not in the operators themselves. 

One then observes the following recursive operator identities: 

crt = \crl[ai,[H,al]], (8) 

< = - <2 - ^/-iM-i]' J = 3, 4 . . . , n, (9) 

which generate the entire set {o"/; j = 1, . . . , n} starting from just H and af. Similarly, 
= l,...,n} are generated by Hermitian adjoints of (!8|9|) . Clearly, the set 
{cr^, <yj',j = 1, . . . , n} then generates, by multiplication and addition, the entire algebra 
B{J-'n)- We note that we can also take only a pair up-flip and down-flip Lindblad 
operators from the opposite edges, say and a~, i.e. if a = 5 = (or, alternatively, if 
/3 = 7 = 0), since the operator recurrence ( P[9i) . or its Hermitian conjugate, can also be 
started at j = n (i.e. it is symmetric with respect to the replacement j ^ n + 1 — j) . 



3. Perturbative (weak coupling) solution 

We start by considering the weak-coupling regime, where all the four rates are small, 

a = ea, (3 = eb, 7 = ec, 6 = ed (10) 

and e is considered as a small parameter. Let us write the NESS density operator as a 
formal power series 



00 



Poo = 5^(ie)V^). (11) 
p=0 

Plugging the ansatz (JTTl) into the flxed point condition (jl]) results in an operator valued 
recurrence relation for the sequence {p^^^;p = 0, 1, 2 . . .} 

^ -I?op(^-^\ ifp = l,2,... 

where Vn = oD- + bV^ + cD- + d'D„+, so that V = eVn. In Ref. il4j we have shown 

u u u u Jl ' c a 

that for a particular case of symmetric driving 

^ , _ 1 -P , ^ _ 1 + P /,oN 

'^sym — '^sym — ^ ; "sym — Cgym — „ ; \ '^l 
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one can express the zeroth, the first, and the second order of the perturbation series 

2>Sin = 1, (14) 

2Vi^L = M^-^^), (15) 

rp^Sr^ = ^{z-z^f-^\z^z\ (16) 

in terms of a non-Hermitian matrix product operator 

Z= J2 (L|A,,A,,---A,jR)a^^®a^^---®a^", (17) 

(si,...,s„)G{+ -,0}" 

where 0"° = 1. Aq, A± is a triple of near-diagonal matrix operators acting on an auxiliary 
Hilbert space Ti spanned by an orthonormal basis {|L), |R), |1), |2), . . .}: 



|L) (L| + |R) (R| + ^ cos (rA) |r) (r | 



r=l 

oo 



A+ = |L)(l|+cf^sin(^2 ^ 



r=l 

oo 



A |r)(r+l|, (18) 



|1)(R| -c-i^sin ((2 [^J+l) a) \r + l){r 

r=l 



where A = arccos A G M U iM and \_x\ is the largest integer not larger than x. Constant 
c G C — {0} is arbitrary, but it is perhaps suitable to choose c = 1 for |A| < 1 (A G M) 
and c = i for | A| > 1 (A G iM) making the matrices f lTSj) always real. The key property of 
the Z operator which is proven and used extensively in [H] is the almost-commutation 
(or conservation law) property 

[H,Z] = -al + a:. (19) 

We shall proceed to show now that this solution can be used also to express 
(perturbatively) the leading orders of NESS density operator for arbitrary drivings 
(arbitrary rates a,b,c,d). This we will do by writing quite general ansatze for the 
zeroth and the first order 

2"p(0) _ (^0 ^ ^^z^0n _ ^20) 

2V^^ = - Z^)R, (21) 

where u (related to an average magnetization in the zeroth order u = {erf) and /i 
are still undetermined parameters. We note that the operator R can be expressed in 
terms of an exponentiated total spin projection. 



R = Vl - z/2exp{(artanhz/)M'}, = ^a^', (22) 

hence it can be shown to commute with the Hamiltonian ([1]) and the Z-operator (ITTj) 

[H,R] = 0, [Z,R] = 0, (23) 

t [Z, AP] = follows from the fact that all the Pauli terms of ([TTl) contain the same number of cr+ and 
(7^ factors. 
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guaranteeing the zerorth order p = condition f ll2p and making the ordering of the 
terms in the first order fl2T]) not important. 

Plugging the ansatze ( 1201I2T]) into the equation [H^ p^^^] = —T)p^^^ and using the 
equations fll9f23p one gets an identitiy jjia\ — fia^ = j^{{b — a — z/(a + b))al + {d — 
c — i/(c + d))a^ — vib — a — i/(a + b) + d — c — ^{c + d))l}, which immediately fixes the 
unknown parameters z/, fi. Namely, the average magnetization reads 

. = ^^4^^, (24) 
a+b+c+d ^ ' 

and the effective driving is 

2 (be — ad) 

^ = (l-u2^(a + b + c + d)- ^^^^ 

We note that the ansatz for the second order p^^^ (fT6!) cannot be extended to general 

asymmetric boundary conditions in a similar way as the zeroth and the first orders 

Jo,i) 

j'sym 



4. Non-perturbative (extreme driving) solution 

Let us now focus on a nonperturbative (exact) solution of the case with just two Lindblad 
channels 



L^ = ^l3at, L2 = ^a- (26) 

For the symmetric situation /3 = 7 =: e this corresponds to the extreme driving case 
p = 1 of ( |T3i) which has been solved in exactly in Ref. |15], namely 

sym _ (27) 

and Sn is a non-Hermitian matrix produet operator 

Sn= (0|A; A;^---A;jO)a^^®a^^---®a^" (28) 

(si,...,s„)e{+,-,o}" 

where (T° = I2 and Aq, A'_j_ is a triple of near-diagonal matrix operators acting on 
an infinite-dimensional auxiliary Hilbert space %' spanned by an orthonormal basis 
{|0),|1),|2),...}: 

00 

A[, =|0)(0| + ^a°|r)(r|, 

r=l 
00 

A; = i£|0)(l| + ^a+|r)(r + l|, (29) 



r=l 

00 



A'_ = |1)(0| + ^a;|r + l)(r 



r=l 
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with matrix elements (writing again A = arccos A) 

/ . sin(rA) 

-= cos (r A) + ie ^ ' 



'Hk-i ~ "^si^ (2fcA) + ie 



2 sin A 

csin((2A;-l)A)sin (2A;A) 



2(cos((2A;-l)A) + r2fc_i)sinA' 



+ • roi \\ ■ c(cos (2fcA) + T2k) , . 

«2fc = csm (2A;A) - le -— , (30) 

/ sm A 

_ _ sin((2A;-l)A) . cos ((2A;-1)A) + rafc-i 

sin((2A; + l)A) . sin (2A;A) sin ((2A; + 1)A) 
c 2c(cos(2A;A) + r2fc) sinA' 

Constant c G C — {0} and signs G {±1} are arbitrary, i.e. all choices of c, give 
identical operator Sn (l28l) . 

The crucial ingredient in the proof of Ref. [15] was the following recursive identity 
satisfied by operators Sn 

[H, Sn] = -ie{a^ ® Sn-i - Sn-i ® a^), (31) 

which can be understood as a 'non-perturbative' analog of the commutator relation (fT9|) . 

We shall now show how one can generalize the ansatz (12711281) in order to incorporate 
the asymmetric driving fl26l) for any rates /3, 7. We start by recognizing the following 
non-unitary symmetry of the XXZ dynamics: 

Lemma: Let v G (—1, 1) he a real parameter and V\ : — a non-unitary, 

but non-degenerate linear map a set 0/ 2 x 2 matrices onto itself, which is completely 
specified by its action on the Pauli basis a'^ = Vi(cr'*), s G {0, +, — ,z}.- 



a" = a° - 



^/(T^ (32) 



1 



= -{a^ - va 



Vi induces a one-to-one linear map V = Vf""^ of n- spin Pauli algebra B{J-'n), which is 
completely specified by V{cr^^ (g) a*^ ■ ■ ■ (g) o"*") = 5"^^ 0"^^ ■ ■ ■ (g) 5"'^". Then, V commutes 
with the Heisenberg dynamics of the XXZ chain (Q]), i.e. for any x G B{J^n)- 

[hMx)] = V{\H,x\). (33) 



Proof: Writing the Hamiltonian ([T]) as a sum of two body terms H = Yl]=i 
a sufficient condition for (155]) to hold is V(a;)] = Due to 

the linearity of this relation in x, the latter can be considered to be of the form 
X = (g) xyj+i] ® X[j+2,n], where the first, or the last, factor are taken as trivial if 

j = 1, or j = n — 1, respectively. Proving the lemma is then equivalent to showing ( l33i) 
for n = 2, i.e. 

[h,a'^a')]=V{[h,a'^a']), s, t G {0, +, -, z} (34) 
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where h = 2(T+ (g) ct + 2cr 0"+ + Aa^ (8) a^. This follows from observing: 

[h, a+ ® a°] = 2A(T+ (g) a" - 2a" ® a+, 

[/i, a+ ® (t"] = 2Aa+ ® (t° - 2(t° ® (t+, 

[/i, cr+ (g) cr"] = 0-° (g) a'' - 0-'' (g) cr°, 

[/l, cr° (g cr"] = 4(7+ (g - 4(7" (g (T+, 

[h,a'^a'] =0, for sG{0,+,-,z}, (35) 

and related Hermitian conjugated identities and identities with swapped tensor factors, 
and checking exactly identical identities with cr* replaced by 5"*. QED 

We note that the non-unitary symmetry map V is not a canonical transformation, 
namely o"** do not satisfy the same commutation relations as the Pauli matrices a*. 

Now, we are in position to state the main result: 

Theorem: A unique NESS of the flow ^ with ( flf^) can be written as 

where R = (a^+ua^)'^"' is the (unnormalized) weak-coupling limit ^2^), and the operators 
Sn are given in terms of the same form of the MPO as in the symmetric case ^2^1 

~Sr.= 5^ (0|A;^A;^-..A;jO)a^^®a^^---®a- (37) 
{si,...,s„)e{+,-,o}" 

but in a different basis [3^] . i.e. Sn = V{Sn), with the magnetization parameter: 

and the effective coupling e (to be used in explicit expressions (dPjE^) for the triple of 
matrices A'_^_, A'_, Aq): 

s - (39) 

/3 + 7 

Proof: Uniqueness of the fixed point poo has been shown, also for the restricted case of 
driving fl26|) , in section |2j 

We need to show that ( 136|1 satisfies the equation (jll), which, together with [_R, S'„] = 
again following from the fact that all the Pauli terms of ( !37|) contain the same number 
of and a~ factors, is equivalent to the condition 

\\HXs\},=V{SnSlR)R-^ (40) 

Note that, since S^ = V{Sn), the Lemma implies 

[H, Sn] = -is ((a" - z/(T°) ® Sn-i - Sn-i ® (a^ - z/a")) , (41) 

where e = e/(l — z/^) and e is still unspecified. Following the idea in the proof of Ref. [T5] . 

we can use the algebraic relations among A'^ to write S'„ =: cr° (g Sn-i + (g P, 



n-l 



Sn-1 ® + Qn-1 ® cr , defining Qn-i and Pn-i as some linear operators over J^n~i- 
This merely expresses the fact that in (5*^ the first always comes to the left of all 
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a , and the last a comes to the right of all cr+. Straightforward calculation, employing 
f HTj) and definitions of Pn,Qn, then expresses the left- hand- side of fHOj) as 

- 25„_i5Li ® - ^^^°) + (^^ - l)Qn-i5l_i ® a- + (z/ - 1)5„_iQLi ® ^+}- (42) 
For the right-hand-side of pUj) . we use the fact that the map 

p V{pR)R-^ (43) 
nontrivially acts only locally, on sites j = 1 and j = n, finding 



- ^^Sn-iSi_^ ® (a" - w°) - 7Qn^i^l-i ® - 7S„_igl_i ® a+. (44) 

Indeed, these two expressions, (H2|) and (l44l) . are linear combinations of exactly the same 
terms. Comparing term-wise results in conditions e = (3/{l + u) = 7/(1 — z/), yielding, 
with e = {1 — i'^)£, exactly the expressions fl38l) and fl39l) . QED 



5. Conclusion 



Computation of expectation values of physical observables in NESS, based on our new 
results for generalized boundary conditions, can be facilitated in terms of the transfer 
matrices which can be constructed exactly as in Refs. [T4l [15]. The effect of the 
asymmetric boundary conditions is essentially an offset of average magnetization, which 
is a consequence of the non-unitary symmetry V described in the Lemma of section |H 
Inspiring discussions with W. de Roeck, H. Spohn and F. Verstraete, as well as 
the hospitality of NORDITA in Stockholm where a part of the results reported here 
were obtained, are warmly acknowledged. The work has been supported by the grants 
Jl-2208 and Pl-0044 of ARRS (Slovenia), and ESF (supporting NORDITA program). 
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